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We consider theoretically tunneling characteristic of a junction between a normal metal and 
a chain of coupled Majorana bound states generated at crossings between topological and non- 
topological superconducting sections, as a result of, for example, disorder in nanowires. While 
an isolated Majorana state supports a resonant Andreev process, yielding a zero bias differential 
conductance peak of height 2e 2 /h, the situation with more coupled Majorana states is distinctively 
different with both zeros and 2e 2 /h peaks in the differential conductance. We derive a general 
expression for the current between a normal metal and a network of coupled Majorana bound states 
and describe the differential conductance spectra for a generic set of situations, including regular, 
disordered, and infinite chains of bound states. 



Topological materials are of large current inter- 
est, in part because of their potential for topological 
quantum computing and their interesting non-Abelian 
quasiparticlesP One variant of this is topological super- 
conductors where the low en ergy quasiparticles in ad- 
dition are Majorana FermionsPEl Currently, there is an 
active search for materials that can host such particles, 
either in certain p-wave superconductors, or semicon- 
ductors with proximity induced superconductivity and 
strong spin-orbit couplingl^E' Because it takes two Ma- 
jorana Fermions to form a usual Dirac Fermion that 
can couple to other degrees of freedom, detection of the 
state of the Majorana fermion system requires non-local 
measurements or interferometryPHfll In contrast, a local 
tunnel current, being independent of the parity of the 
topological superconductor, does not reveal information 
about the state of the Majorana fermions. Nevertheless, 
a tunneling probe coul d detec t the presence of a Majo- 
rana bound state (MBSpAM^and the detection of Majo- 
rana bound states is the first major challenge in this field. 
Tunneling contact to an isolated Majorana state give rise 
to a resonant Andreev process that gives a zero bias con- 
ductance peak of 2e 2 //iP^ With two coupled Majorana 
states cross correlations of the current into eac h cou ld 
also show their existence and non-local character! 13 ! 14 ! 

However, because of material difficulties, isolated Ma- 
jorana states might be rare. Rather it is to be expected 
that density fluctuations will generate a random con- 
figuration of topological/non-topological boundaries, at 
which Majorana states will be located. For strong dis- 
order the distances between these states are sufficiently 
short for the MBS to overlap and therefore it is impor- 
tant to understand how a network of coupled Majorana 
fermions maps onto the tunneling characteristic. This 
problem was recently considered in Ref. [TS] in the weak 
coupling regime, using a renormalization group to reduce 
a chain to a sum of single Majorana pairs on a logarith- 
mic energy scale. 

In this paper, a theory for tunneling between a metal- 
lic probe and a collection of coupled Majorana states in 
the strong coupling regime is developed, and experimen- 
tally relevant situations are addressed. This is done in 
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FIG. 1. (Color online) Illustration of a semiconductor with 
induced superconducting order parameter from an adjacent 
superconductor (S). Spatial variations in density or supercon- 
ducting order parameter create crossings between topological 
(TS) and non-topological segments when the density crosses 
the critical density (vertical (red) lines). Majorana bound 
states (MBS) (circles) are located at each crossing point and 
the distances between them determine the coupling matrix 
elements iy of the resulting Majorana network. A tunnel- 
ing contact (N) probes the network by tunneling into the end 
Majorana mode, with tunneling density of states V. 

the limit where the voltage eV, the tunneling broaden- 
ings r, and the hopping matrix elements between any 
two MBS all are much smaller than the superconducting 
energy gap A. The regime of large A is well suited for 
characterization and detection of the MBS, because in 
absence of Majorana states the Andreev conductance is 
of the ordei!^ (e 2 /h)(T / A) 2 and thus much smaller than 
the resonant Andreev current carried by the Majorana 
states. Examples for different number and configurations 
of Majorana states is given and the case of a uniform infi- 
nite chain is solved exactly. Finally, disordered Majorana 
chains arc addressed. Disorder is introduced as random 
nearest neighbor couplings and it is show to reduce to a 
finite chain, truncated by the first weak link (quantified 
below) in the chain. 

The borders of the topological superconductor seg- 
ments give rise Majorana bound states. These states are 
zero energy solutions to the Boguliubov-de Gennes equa- 
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tions for the geometry in question. The general form of 
a MBS is 



It = 



(i) 



The Majorana Fermion has the properties that r y i = 7] 
and 7j 2 = 1. The superconductor Hamiltonian, describing 
the coupled Majorana state network, is 



(2) 



The tunnel Hamiltonian between the normal metal and 
the superconductor is 

ka 

where Ck,a are lead-electron annihilation operators and 
^a{ x ) the superconductor electron- field operator. As ex- 
plained above, for (eV, r)«A the Majorana states con- 
tribution to the current dominates. Using Nambu rep- 
resentation, \& = ("^f, ^j,, "^j, the projection of the 
field operator *f? onto the manifold of Majorana states is 

« Ei7i(/ t) i(*).4i(a:) ./iiW./f.iW)' which then 
leads to the effective tunnel Hamiltonian describing the 
coupling between the lead and the Majorana states 



(3) 



ka i 



where Vfe CT .i = J dx f a ^xjt^fx). The current operator is 
given by the rate of change of the number of electrons in 
the normal lead 

I = = ~ie[H Ti N]/h = ^ Re]T (v^Gf^O) 

kai 

(4) 

where the lesser Green's function combining ka and i is 
defined as Gf fe(T 



Gfka 



(<) = i ^Cfe CT 7i(i)^) which is written as 

(') = £[■ 



(5) 



where Gy(r,r') = -i (T^r^r')) is the full Keldysh 
time-ordered Green's functions for the Majorana opera- 

tors, and = (r(c fcff (r)cL(r')) o is the un- 

perturbed normal lead Green's function. By choosing 
the chemical potential of the superconductor as a refer- 
ence, the general current formula is derived to be (see 
Appendix) 

I = t J dej M(oj) [f(-oj + eV) - f(u - eV)] , (6) 

with / being the Fermi-Dirac distribution and 

M(lu) = Tr [G r {oj)T* (-u)G A (u)T{u)] . (7) 



Here the retarded Majorana Green's function is 

G* = 2(w-2it + i(T u + r*_ u ) - (A^-A^)) -1 , 

.. (8) 

where t is an antisymmetric matrix, while the Hermitian 
matrices T and A are 

r tj (w) - 2tt Vk^J (w - e ka ) , (9) 



ka 



duo' Tij(uj') 
2tt u) — u)' 



(10) 



If the coupling matrix respects particle-hole symmetry, 
r(w) = r*(— uj), the current is antisymmetric I(V) = 
— /(— V) (see Appendix for more details). 

The expression ^ is a general finite temperature ex- 
pression for the current into a Majorana state network in 
terms of matrices describing the coupling to the normal 
lead and the Majorana network. The general formula is 
straightforwardly extended to the case with more normal 
metal contact connected to the network.^ 

If the Majorana bound states are separated in space 
by a distance much longer than the normal metal Fermi 
wavelength, off-diagonal terms of E^ will tend to average 
out due to the fast variation of the phase of V ka i . In this 
case, it is a good approximation to set r^(w) ~ SijTa(uj), 
which is assumed from here on. Moreover, assuming a 
weak energy dependence Vkai arL d hence a constant Ta, 
so that Ay = 0, (so-called wide-band limit), the differen- 
tial conductance reduces to 



dl _ 2e 2 
dV ~ ~h 

with 



df(u - eV) 
du> 



G R (co) 



2it + i2T] 



(11) 



(12) 



and r is a diagonal matrix. 

For a single isolated Majorana state with tunnel broad- 
ening coupling the Green's function is: G R = 2/(w-M2r) 
and the zero temperature differential conductance is eas- 
ily obtained as 



dl 2e 2 



dV 



4r 2 



(13) 



which confirms that the resonant Andreev tunneling with 
zero bias conductance G = 2e 2 /h^ With two Majoranas 
coupled by tunneling t and only one of them coupled to 
the lead the differential conductance is 



dl 
dV 



2e^ 
IT 



(2eVr) 2 
((eV) 2 - 4i 2 ) 2 + (2evry 



(14) 



which has a dip at zero voltage and peaks at eV = ±2t, 
where the conductance again reaches 2e 2 /h. In fact, a 
very general statement holds for tunneling into the end 
of a chain, namely that with an odd number of coupled 
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FIG. 2. (Color online) The conductance for tunneling into 
the end of a chain with 1,2,5, and 8 Majorana states coupled 
by t — r. A general feature, valid also for disordered arrays, 
is that for even number of coupled Majorana modes, the con- 
ductance is zero at zero voltage, whereas for an odd number 
it is given by 2e 2 jh. Moreover, the conductance has in general 
n— 1 zeros, where n is the number of Majoranas in the cluster 
that couples to the probe. The dashed (red) curves show the 
result when the two first Majorana states both coupled to the 
lead, the second one with strength F22 = fn/2. 



Majorana states the zero bias conductance is 2e 2 jh, and 
with an even number the zero bias conductance is zero. 
This can be shown by the inversion in Eq. (12 1 setting 
uj = and for an arbitrary chain matrix tij. 



Moreover, 

for a cluster with n MBS the differential conductance ver- 
sus bias voltage has n — 1 zeros and n voltages where 
dl/dV = 2e 2 /h. If the normal metal electrode overlaps 
with more than one MBS these conclusions change, as 
shown in Fig. [2j where the conductance for some exam- 
ples is plotted. 

With many coupled MBS in the chain the conductance 
oscillates between 2e 2 /h and 0, as seen in Fig. [2] As the 
number of sites in the chain is increased, the period of 
the oscillations decreases. If the period is smaller than 
temperature the conductance will average to a value be- 
tween the two extremes. The same occurs for an infinite 
homogeneous chain, which is considered next. 

With an infinite chain of Majorana states with nearest 
neighbor couplings 
MBS is Gfi = 2g u , where 



t^ , the Green's function for the first 



5n 



1 

-4|ti 2 | 2 5 22 ' 



(15) 



and where g 22 is the Green's function for the network 
starting with site 2, decoupled from site 1, and where 
(ffii) -1 = w + 2«r. An illustrative example is a homo- 
geneous chain, i.e. with all couplings identical tij = t. 
Then the Dyson equation for g 22 is 



1 



it 2 



922 



uj + in uj + in 



533522- 



(16) 
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FIG. 3. (Color online) The conductance for tunneling into 
the end of an infinite Majorana chain with identical tunneling 
couplings. The parameters \t\/T ranges from 5 to 1/8, from 
outside in. 



Since all connections are equal g 22 = <? 33 , and hence 
{4:t 2 /uj)g 22 — g 22 + 1/uj = 0, which gives can be solved 
for g 22 . Choosing the correct branch cuts^ and setting 
this into the Green's function (15), the differential con- 
ductance is derived to be 



dl _ 2e 2 
dV ~ IT 



4r (4r+ % /(4t) 2 -( e y) 2 ) 

( e y) 2 +(4r +v /(4t) 2 -( e v) 2 )^' 
(4r) 2 

(\eV\ + ^J(eV) 2 -(4t)A 2 + (AVy ' 



\eV\ <4|t|, 
\eV\ >4|t|. 



(17) 

This is an interesting expression with a line shape that 
strongly depends on the ratio t/T, which is shown in 
Fig. [3] Furthermore, for eV = it reduces to 



dl_ 
dV 



2e 2 



2T 



v=o 



h t\ + 2^ , 



(18) 



which shows that for small tunnel broadening compared 
to the bandwidth of the chain the zero bias conductance 
is 2e 2 /h, which was expected because it corresponds to 
tunneling into an effectively isolated Majorana state. 

As the last situation, which might also be the most 
experimentally relevant, we now discuss different realiza- 
tions of long disordered chains, sampled for example by 
scanning the average density and creating a different set 
of crossings of the topological superconductor thresholds, 
as illustrated in Fig. [TJ One could in this way study the 
average conductance of a network, given a distribution of 
nearest neighbor MBS couplings Uj. The tunneling cou- 
pling between two neighboring MBS depends both on the 
distance between them and the deviation from the crit- 
ical value for the topological/non-topological transition, 
with exponential dependence on both, as was shown by 
Shivamoggi et alP^ for a specific example. The distribu- 
tion of tunneling couplings is thus a complicated convo- 
lution of amplitude fluctuations and the level-crossings 
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FIG. 4. (Color online) Both panels show the conductance for 
a chain with 10 sites and one weak link with tunnel coupling 
iwcak = 0.25r, while the other links have t = T. The weak 
links are the 1-2 and 5-6 connections for the left/right panel, 
respectively. The thick (blue) curve is the conductance for 
kT = O.ir, the thin gray curve for kT = 0, while the crosses 
show the conductance for the chain truncate at the weak link 
also with kT — 0.1F. It is clearly seen that the truncated 
approximation works well, because the chain after the trun- 
cations leads to structure barely resolvable, because fesT is 
not much smaller than the width (ss 0.125) of the additional 
resonances. 



statistics,^ in this case with two crossings. The result- 
ing distribution of tunneling couplings is an interesting 
problem in itself. However, instead of pursuing this line, 
we focus at the generic behavior expected for a given 
configurations of the tunneling couplings in the chain. 

As we learned for the infinite chain, different behav- 
iors occur depending on the ratio of T to the tunneling 
couplings tij. For the random chain the ratio of T to 
the spread of tunneling couplings turns out to be cru- 
cial. Clearly, if the spread in tunneling couplings is much 
smaller than F, the average conductance resembles that 
of the homogenous infinite chain, which we have verified 
by numerical simulation. 

In contrast, with large fluctuations in the tunneling 
couplings, the infinite chain will be effectively truncated 



into a finite chain, where one of the tunneling coupling 
happens to be much smaller than ksT. To see this, invert 
the matrix in Eq. ( 12 ) and pull out the dependence on 



the weak link and write is as 



G R {ll,w) 



2Do 



-Dl,oo — 4t n n+1 -Dl.n-Dn+l.oo + iTD2,oo 

where i wea k = in,n+i is the weak link, and where Dij is 
the determinant of the matrix (w — 2t) for the isolated 
chain between sites i and j, but with t n ^ n +\ = 0. The 
weak link has two effects, 1) it gives small shifts of the 
existing resonances and 2) it creates new resonances. The 
new resonances, however, have widths that scale with the 
square of the weak coupling r wca k cx i^, oak / (^)r, where (t) 
denotes typical couplings in the first part of the chain. 
Therefore the new resonances introduced by the chain 
after the weak link is not resolved if kst > r wea k. An 
example of this is shown in Fig. [4j where the conductance 
of a truncated chain is compared with that of a full chain. 

In summary, the differential conductance for a junction 
between a normal metal and topological superconductor 
hosting a network of Majorana bound states has been 
studied. Different configurations of the interacting net- 
work of bound states give rise to distinct tunneling spec- 
tra. Long chains with fluctuating tunneling couplings 
is truncated into a finite chain once a coupling becomes 
smaller than a certain critical value, determined by tem- 
perature. 
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Appendix A: Current formula for a Majorana state 
network coupled to normal-metal electrodes 

The low energy model Hamiltonian projected onto the Ma- 
jorana subspace is 



The corresponding Keldysh contour Green's functions are 



H — Hn + Hs + Ht 

hn = y ] £ p cpc p , 

ka 

H s = |E^ 7i7 i' 

ij 

H T = ^(Vpicl - V pi c p )~/i, 



(Ala) 
(Alb) 

(Ale) 

(Aid) 



where p = k,a. In the non-equilibrium case the voltage on 
the normal metal is eV. The coupling matrix tij is real and 
obeys tij = — The current operator is 



/ = -eN = -ie[H T ,N]/h 

= Jl E { V pi C lli - VpiH c p) 



= \ E (KGf P (0, 0) - V pi G<(0, 0)) , (A2) 



and 



G<(t,t') = i(ct(t'hi(t)), 
G<{t,t') = i{ li {t)c p {t')). 



(A3a) 
(A3b) 



G ip (T,T') = -i(T K ( li (T)cl(r'))), 

G pi {T,T') = -i(T K {c p {T)>yi(T'))). 

Below we will also need 

G pi (T ) T') = -i(T Jf (4(T)7i(T'))), 

the Majorana Green's function 

Gi j{T,T') = -i <T K ( 7l (r) 7j (r'))), 
and the free lead electron Green's functions 
G p(r,r')^-i{T K (c p (r)cl(r'))) () , 
G p °\ry) = -i{T K (cl(r)c p (r'))) o , 

where the unperturbed expectation value (-} is 
situation with Ht = 0. The self energies 

^y) = E^ G p(^')> 

4(^) = EW^ 0) (^'), 
p 

will also appear. 



(A4a) 
(A4b) 



(A4c) 

(A4d) 
(A4e) 

for the 

(A5a) 
(A5b) 



Now return to the mixed Green's functions G ip , G pi , and G pi . From diagrammatics, direct expansion, or equation of motion, 
we obtain 

G ip (r, t')V£ = E / dT " G *i( r > t ")Z%(t", t') = J dr" Tr [G(r, r")S e (r", r')] , (A5ca) 

p'i ij 

Y,V pi G pi (T,r') = J2j dT " ^tj(r,r")Gji(r",T') = J dr" Tr [S e (r, r")G(r", r')] , (A5cb) 

p'i ij 

E Vp\G pi {T, t') = - E / dr" S«(t, r")G 31 (r", t') = - J dr" Tr [e"(t, r")G(r", r')] ■ (A5cc) 



With this the current becomes 

Next, let us find the Dyson equation for the Majorana 
Green's function. Its equation of motion is 

^G !3 (t,t') = 25 i jS(T,r') + i(T K ([H, Ji ](Thj(T'))) . (A4) 



The factor of 2 because {7i,7 3 } = 2(5; 3 . The commutators: 



[Hom] = ^ E^"[7i'7j',7i] = *E ~ ^)7? = -2iE^'7j' 



i'j' 



[-ffr,7»] = E [(^4 - ^pj c p)7j>7*] = 2E( V ^4 - ^ c p)- 



(A5) 
(A6) 
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Again the factors of 2 come from the unusual commutation relation. With this the equation of motion becomes 



id T G l3 {ry) = 25 i j6(T,T') + 2iY^t ij ,G r j(T,T') + 2^2(V pi G pj (T,T') - V£G pJ (t, t')) . 

f P 



By Eq. ( A5c I the Dyson equation closes 



ia T G i3 (r,r') = 2<5y5(r,r') + 2i V Gi'i (t, r') + 2^/ dr " ( E ^'( r > r ") + E ^(r,r")) G j0 (t 



(A7) 



(A8) 



or in shorthand notation 

(id T - 2it - 2E) G = 2, 
which has the solution 



G = G° + G°EG, 



(A9) 



(A10) 



where the unperturbed Majorana Green's function is 

(id r - 2it) G° = 2, (All) 

and where the Majorana self energy is 

E = E e + S' 1 . (A12) 

The same result can be derived using diagrammatics, in 
which case the factors of 2 then comes from the 2 ways 4 
Majoranas can pair. 

The retarded components of the self energy are 

^(«)=Ei^=^+M«), (A13) 
v 

^M=E^=^-M-), (AM) 
p 

and as usual S A = Here the T and A matrices are 

defined as 

r M (w) = 2tt £ V^V^tf ( W - e„) , (A15) 

V 

and they are both Hermitian matrices, so that Y1 R — S A = 

i(r w + r* w ) . 

The lesser components are 



(A17) 



E^(w) = iT£ (-w) (1 - /(w - el/)). (A18) 



Now go back to current and use that (_BG) < = B R C < + 
B < C A to get 



e f dui 



Ti- 



ft./ 2tt 
In equilibrium G^' eq 



[(g5-g^s;, -i-g. (s:.: 1 -e' r 



(A19) 

■i (G„ — G„) / w and the current 
is zero. The Majorana lesser function is 



G^J — G^E^G^ — iG R (Tufu-ev + r_ w /- 
since f^-ev = 1 — f-u+eV- Together with 



ui + eV 



(A20) 



G£ - GZ = G* (£* - E^j G^ = »G£ (T„ + r' „) G^ 

(A21) 

this finally leads to 



I=- dw M(w)(f-^ +eV - U- eV ), (A22) 



with 



M(w) = Tr [G^T* 



iG,,, r, 



w u; w j 



(A23) 



which is the final general result for the current from a nor- 
mal metal into a network of Majorana fermions. 
The retarded Green's function is 

G* = 2 (w - 2it + i (T u + r* u ) - (A w - A* J)" 1 . (A24) 

In the electron-hole symmetric case T u = r* u = r^„, 
and hence A w = — Al w , it follows that 



G" =- G" = — ( G 



(A25) 



and therefore 
M(-u) = Tr [g^KG^T 

= Tr [(r^G^rl„G^) T ] = M{u) 



= Tr |G^ T riI,G^ T r 



1} 

(A26) 



Therefore, if the leads do not break electron-hole symmetry 
the current is anti-symmetric I(V) = —I(—V). 



